Coexistence of long-range and algebraic correlations for short-range valence-bond wave functions 

in three dimensions 
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We investigate nearest-neighbor valence-bond wave functions on bipartite three-dimensional lattices. By 
performing large-scale Monte Carlo simulations, we find that long-range magnetic order coexists with dipolar 
four-spin correlations on the cubic lattice, this latter feature being reminiscent of the Coulomb phase for classical 
dimers on the same geometry. Similar properties are found for the lower-coordination diamond lattice. While 
this suggests that the coexistence of magnetic order and dipolar four-spin correlations is generic for such states 
on bipartite three-dimensional lattices, we show that simple generalizations of these wave functions can encode 
different ordering behaviors. 
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Given the difficulty to obtain ground states for even sim- 
ple model Hamiltonians, one approach in trying to find novel 
strongly correlated states is to study variational or trial wave 
functions directly, and to seek the Hamiltonians which they 
minimize in a subsequent step. Most prominently, such an ap- 
proach has been successfully pursued in the case of Laughlin's 
wave function for the fractional quantum Hall effect yj] . 

For unusual magnetic states, wave functions defined in 
terms of valence-bonds (VBs) have an illustrious history 
J2, Etl- In their simplest incarnation, in the context of quan- 
tum dimer models, a dimer represents a VB. This has estab- 
lished a connection between the venerable subject of close- 
packed hard-core dimer coverings fl-Q] and the physics of 
strongly correlated electrons, and has led to important con- 
ceptual progresses, including the discovery of a liquid phase 
with exponentially decaying correlations ]8|]. More recently, 
the problem of dimers on the cubic and other bipartite three- 
dimensional (3d) lattices has received considerable attention 
JsKIol]. The hard-core constraint for dimers allows for the for- 
mulation of an effective theory in terms of a coarse-grained, 
divergence-less, magnetic field fl, implying that dimers on 
bipartite 3d lattices realize a Coulomb phase Hill , charac- 
terized by dipolar dimer-dimer correlations. We also briefly 
mention the close relationship [Till between hard-core 
dimers on the diamond lattice and spin-ice systems on the 
pyrochlore lattice II 1 311 - that have been the focus of intensive 
research in recent years I114[ Il5 1 . 

VBs or related degrees of freedom also appear in several 
other fields of research, including quantum chemistry, from 
where they originated (see e.g. Refs. [l6| and 1 1 7b . The inti- 
mate relationship between superimposed VB configurations 
and loop models was later established, the latter being inten- 
sively studied in statistical mechanics 11811 . for e.g. their geo- 
metric transitions II 19^2 ill . In quantum information, entangle- 
ment properties of short-range VB wave functions have been 
recently studied I22l - l24ll . Finally, we mention the recent ex- 
perimental realization of a plaquette VB state with ultracold 
atoms loaded into optical superlattices [25]. 



Within the framework of quantum magnetism, we remark 
that hard-core dimer approaches, while leading to exact re- 
sults even in the 3d case, do not a priori incorporate actual 
spin degrees of freedom. However, it is also crucial to be able 
to account for magnetic order and its interplay with phases 
lacking it, in order to obtain a full picture of the nature of 
possible states. In this work, we analyze the effects due to 
the incorporation of spin degrees of freedom by characteriz- 
ing nearest-neighbor VB (NN-VB) wave functions defined on 
the 3d bipartite simple cubic and diamond lattices (see Fig . [TJ . 
Such states can be regarded as 3d generalizations of the NN- 
VB wave function on the square lattice originally introduced 
by Sutherland 12611 . that has been recently studied via Monte 
Carlo (MC) simulations 123, H- NN-VB wave functions on 
the cubic lattice have been previously analyzed in Ref. |29|], a 
combined mean-field and MC study of a family of Hamilto- 
nians with interactions of varying ranges for spins of varying 
lengths. This work already found the presence of magnetic 
order - see below (we also mention recent studies of NN-VB 
states on frustrated 2d lattices |24 1|]). Here we present, 
for cubic and diamond lattices, a thorough evaluation of the 
correlations - e.g. the presence of dipolar four-spin correla- 
tors - encoded by NN-VB states, characterize the ensemble of 
loops encountered, and provide an outline of a field-theoretic 
description of an extended family of NN-VB wavefunctions. 

We consider the NN-VB state defined as an equal- 
amplitude superposition of all hard-core dimer coverings c on 
the considered lattice, 



W = £|c>, 



(i) 



where a singlet state, 



f(|t4i) - \Utj)), is attached 



to each dimer joining site i on sublattice A to one of its NN 
sites j on sublattice B. This fixes the overall sign of the wave 
function, so that stochastic methods are applicable J2l- The 
intimate relationship between classical dimers and NN-VB 
states in this definition suggests that the Coulomb phase found 
for dimers may also be encoded in 5f/(2)-invariant NN- 
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FIG. 1. (Color online) (a) Cubic and (b) diamond lattices consid- 
ered in the present work. Thick light-shaded links in (a) represent 
singlets in a NN-VB state [Eq. (T)]. We also indicate special direc- 
tions (with respect to a reference dimer depicted as a thick dark line), 
along which four-point correlations are computed. 



VB wave-functions. It is our main goal here to test this con- 
jecture, by performing state-of-the-art MC simulations. We 
find an intricate interplay between magnetic and VB (energy) 
correlations, most remarkably in the form of Neel order coex- 
isting with dipolar four-spin correlations. 

Numerical method and observables — A crucial aspect 
of the MC simulations is to take into account the non- 



orthogonality of VB configurations Il26ll . Following Ref . 0, 
we write the expectation value of an observable O in the wave 
function Eq. (Q~|i as 



(ci\0\c 2 ) 



Z ct?i (Cl|C2) 



ci|c 2 ) 



(2) 



Z = stands for the normalization and the overlap 

between NN-VB configurations \ci) and c 2 ) is given by 
(ci|c 2 ) = 2 n £ ( - 1 > 2 )-tt, with N c (l,2) denoting the number 
of loops in the transition graph for both configurations and N 
the number of sites Il26ll . Alternatively, the normalization can 
be written in a pure loop representation as y26 1 



2(C) 



(3) 



where the sum is over the set of all loop coverings obtained 
by superposing two NN-VB configurations, and ni(C) the 
number of loops of size I in the loop covering C. The es- 
timator (ci\0\c2) / (ci\C2) for O in Eq. (0, determined by 
the loop structure in the transition graphs (for instance, the 
squared staggered magnetization m| relates to the squared 
loop-lengths; see Ref.l32l). can thus be sampled by generating 
pairs of dimer configurations, by relying on an efficient worm 
MC algorithm i33ll . Following this procedure, we probe two- 
and four-spin correlation functions for NN-VB states on the 
simple cubic and diamond lattices, for system sizes of up to 
L = 64 (cubic) and L — 32 (diamond, 2 x 32 3 sites) with 
periodic boundary conditions. Note that we do not take the 
existence of topological sectors |@] e xplicitly into account. In 
contrast with the square lattice case 11271 12811 . winding worms 
occur with non-negligible frequency during the simulations, 
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FIG. 2. (Color online) (a) Two-spin staggered correlations as a 
function of distance r for the NN-VB state on the cubic and dia- 
mond lattices (linear size is L — 16). (b) Squared staggered mag- 
netization density versus inverse system size L~ , Extrapolated val- 
ues, obtained from linear fits to the three left-most data points, are: 
(ml)oo = 0.02306(8) (cubic) and (m|)oo = 0.0076(2) (diamond). 



and consequently different topological sectors are sampled. 
This property stems from the existence of long-range antifer- 
romagnetic (AF) order, as we discuss next. 

Magnetic order — We first consider staggered two-spin 
correlations, (— l) r (So -S r ), as a function of distance r for cu- 
bic and diamond lattices with L = 16. As shown in Fig.|2ja), 
correlations seemingly saturate at large distances, suggesting 
that NN-VB wave functions on both lattices display long- 
range Neel order. Previously obtained in Ref. 29 for the case 
of the cubic lattice (see also a related result for higher spins 
13411). this under-appreciated result is rather surprising if we 
recall that such wave functions were originally introduced as 
variational states for magnetically disordered systems. We 
put this observation onto firmer grounds by analyzing the de- 
pendence on system size of the squared staggered moment 

arguments lf35ll predict (m s 
states, in good agreement with our plot Fig. |2jb). Data extrap- 
olation confirms magnetic ordering, with extrapolated values 



L Si j (~ l) nj (^i • Sj). Standard finite-size scaling 
' x to scale as L^ 1 for AF ordered 



(m^oo = 0.02306(8) (cubic) and (m§}oo = 0.0076(2) (di- 
amond). While the lower coordination of the diamond lattice 
weakens AF order, it is not enough to destabilize it. 

Further insight into the occurrence of AF order is given by 
the analysis of the length I of loops appearing in the transi- 
tion graphs (ci|c2). Indeed, the presence of long-range AF 
correlations requires long, system-spanning, loops with a size 
distribution that falls "slowly enough" |32J]. We keep record of 
loop sizes during the simulations, and plot in Fig.[3]the prob- 
ability P(l) for the appearance of a loop of length I for the 
state Eq. (Q]) on the cubic lattice. Intriguingly, P(l) displays 
a two-regime behavior remarkably similar to what was found 
for 'Dirac strings' - the flux loops of the emergent gauge field 
of the Coulomb phase - in spin-ice 112 ill : P(l) ~ l~ a for 
short and P(l) ~ l~P for long loops, with best-fit exponents 
a = 2.51(2) and (3 = 1.07(4). We note that these are the 
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FIG. 3. (Color online) Probability P(l) of observing a loop of length 
I in the transition graph between two NN-VB states when sampling 
the function Eq. (TJ on the cubic lattice, for various system sizes 
L. Two regimes can be distinguished for short and long loops with 
power-law decays with, respectively, best-fit exponents a = 2.51(2) 
and p = 1.07(4) (lines), in agreement 13^1 with the scaling found in 
Ref. I2H1 for a spin ice system. Error bars are smaller than symbols. 



exponents characteristic of random walks 112 111 . The naive ex- 
pectation for the statistics of the nontrivial loops in the tran- 
sition graph generated by the superposition of different VB 
coverings would be that of a self-avoiding but otherwise ran- 
dom walk. Random due to the fact that correlations between 
VBs in the same covering decay quickly (like r -3 ), and self- 
avoiding because the loops cannot self-intersect. However, as 
it is also the case for dense polymer solutions 0711 . the need 
to avoid all the other loops in the transition graph counterbal- 
ances self-avoidance, with the net effect of producing expo- 
nents characteristic of the Brownian random walk. 

We finally remark that the long-range tail of P(l), slowly 
falling as ~ ensures long-range AF correlations 0211 and 
is caused by winding loops in the transition graphs. 

Four-Spin Correlations — We proceed now to the charac- 
terization of 'dimer-dimer" four-spin connected correlations, 



ijkl 



((s < .s J -)(s Jfc -s,))-<s < .s J -)(s fc -s,) 



(4) 



where each pair of indexes i, j and k, I denote NN sites 
(dimers). Anticipating the dipolar structure of these corre- 
lators, we evaluate r 3 aCijki 1211, where r is the distance be- 
tween dimers and a — +1 (— 1) if the sites i and k lie on the 
same (different) sublattice(s), along selected directions high- 
lighted in Fig.Q] 

For the NN-VB wave function on the cubic lattice, we plot 
r 3 aCijki versus r for L = 64 in Fig. Ufa), along the same 
directions [parallel to the vectors e a= o ; i i 2,3 in Fig.QJa)] con- 
sidered in the analysis of classical dimers in Ref.|9l Curves for 
all directions tend to saturate, evidencing a dipolar structure 
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FIG. 4. (Color online) Four-point connected correlations [Eq. {4]l] 
for the NN-VB state on the cubic lattice, multiplied by the sublattice 
sign factor a and the cube of the distance between dimers, r 3 , for: 
(a) L — 64 and specified directions and (b) various system sizes and 
fixed direction eo [special directions are highlighted in Fig.fJJa)]. In 
both panels, symbols indicate data from MC simulations, and lines 
correspond to the finite size expression Eq. ((6). 



where K is a constant and a is the angle between the ref- 
erence dimer along eo and e a= o,i 2,3. Finite-size effects are 
accounted for by the lattice sum expression derived in Ref. S 



aC, 



ijkl 



KL 3 



cosfc y — cosfc 2 ^ ik . r 



(6) 



(5) 



where ki =x ^ y ^ z are the reciprocal vectors of the simple cubic 
lattice. Finite size results obtained from Eq. ©, with the sin- 
gle adjustable parameter set to K = 13.0(5), are depicted as 
solid lines in Fig. [4] and are m excellent agreement with the 
MC data (symbols) for all system sizes and directions, further 
establishing the dipolar structure of four-spin correlations for 
Eq. (G3 on the cubic lattice. 

We have also computed r 3 crCijki for Eq. ([T} on the dia- 
mond lattice, along the directions e a (a = 0,1,2) indicated 
in Fig. QIb)- Fig. |5ja) shows that, similarly to what happens 
in the case of the cubic lattice, four-point correlations are also 
consistent with a dipolar structure [Eq. (0] for the diamond 
lattice, although finite size effects are stronger and satura- 
tion at large r less evident; specially along eo, as shown in 
Fig. f5Jb). We have also performed MC simulations for clas- 
sical dimers on the same geometry, using the algorithm of 
Ref. [3H and find an essentially identical, and equally strong, 
finite size dependence for dimer-dimer correlations also in this 
case [filled symbols in Fig. [2a)]. We have not been able to de- 
rive a finite-size expression similar to Eq. © for the diamond 
lattice, even though we suspect this is possible. 

Discussion — In conclusion, large-scale MC simulations 
of NN-VB wave functions on the cubic and diamond lat- 
tices have allowed us to obtain (numerically) exact results 
for such quantum states, a rare situation for quantum mag- 
netism in 3d. Remarkably, despite of the fact that NN-VB 
states were originally introduced to mimic spin liquid states 
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FIG. 5. (Color online) MC data for r 3 aC ijk i [Eq. ©] for the NN- 
VB state on the diamond lattice as a function of r along special di- 
rections indicated in Fig.[TJb). (a) Data for L = 24 and specified 
directions for classical dimers (filled symbols) and the NN-VB wave 
function (open symbols; data multiplied by 2). (b) Correlations for 
the NN-VB wave function along the direction eo [see Fig.fTJb)] for 
various system sizes. Lines are guides to the eye. 



in 2d, their versions on bipartite 3d lattices are magnetically 
ordered. As a consequence, variational approaches relying on 
NN-VB states, that are only able to describe magnetically dis- 
ordered states in 2d B9I441I1 . may be expected to be useful in 
analyzing the interplay between disordered and SU(2)-broken 
phases on bipartite 3d lattices. 

Additionally, we have shown that four-spin correlations 
possess a dipolar structure, a feature reminiscent of the 
Coulomb phase of classical dimers [9]. This suggests that 
the wave functions investigated here may be employed in a 
variational description of spin-half t/(l) quantum spin liq- 
uids lioll . However, such an approach is hindered by the 
presence of magnetic order, that prevails over dipolar corre- 
lations. An alternative route involves considering VBs with 
enhanced symmetry: for instance, SU(iV) VBs configurations 
have reduced overlaps as N increases. Correspondingly, the 
transition graph loops (c\ \c2) exhibit an increase in the (bare) 
loop tension, leading to a geometric transition where the long 
loops vanish together with magnetic order. Indeed, the large 
N limit corresponds to orthogonal dimer coverings, for which 
spin correlations vanish beyond two lattice spacings. 

Finally, we comment on the basic structure of a field- 
theoretical description of our results, that must account for the 
presence of two distinct emergent Coulomb structures. On 
one hand, four-point dipolar correlations can be expected to 
be described by a fictitious magnetic field, similar to the case 
of hard-core dimer configurations H- In 2d, this approach 
is formulated in terms of a "dimer-height" (generalizing, in 
d = 3, to a "dimer-gauge") variable. Such a dimer frame- 
work has been put in place by Damle et al. 114211 . who pro- 



vide an approximate mapping of the NN-VB problem to a 
classical dimer model with interactions - the possibility of a 
dimer phenomenology having been noted in Ref. |28| (see also 
Ref. |43I) . Combined with the known phase diagram 04411 of 



the corresponding interacting dimer model on the cubic lat- 
tice, the system is predicted to stay in the Coulomb phase 
with dipolar dimer correlations [Eq. (0] as observed. On the 
other hand, this cannot also account for AF long-range or- 
der, as this order can disappear independently of the presence 
of the algebraic dimer correlations. One therefore requires a 
second set of "loop-gauge" variables, whose Coulomb phase 
captures the loop size distribution of Fig. [3] We hope that 
these remarks, together with our numerical results, may stim- 
ulate field-theory work to describe the unusual properties of 
the NN-VB wavefunctions reported herein. 
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